Highly efficient second-harmonic generation from indefinite epsilon-near-zero slabs of 

subwavelength thickness 
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We theoretically predict efficient optical second- harmonic generation (SHG) from a few hundred 
nanometer thick slab consisting of a quadratic nonlinear anisotropic medium whose linear principal 
permittivities are, at the fundamental wavelength, very small and have different signs (indefinite 
medium). We show that, by illuminating the slab with a p-polarized fundamental wave (with 
intensity of a few MW/cm^), a highly efficient scattering of the second-harmonic field occurs when 
the conditions of linear complete slab transparency for the fundamental wave are met. The high 
efficiency of the SHG process, stems from the large non-plasmonic enhancement of the longitudinal 
field, perpendicular to the slab surface, produced by the very small value of the slab dielectric 
permittivities. A suitable nano-structured composite is proposed and numerically designed for 
observing the novel non-phase-matched and highly efficient SHG process from nano-structures. 
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Enhancing optical second-harmonic generation (SHG) 
P, Q is at present one of the most relevant task of non- 
linear optics due to the major role played by frequency- 
doubling in coherent green and blue light sources design 
chemistry [1], biosensing [5|, etc. In situations where 
standard phase-matching or quasi-phase-matching tech- 
niques cannot be used, efficient SHG is generally achieved 
by resorting to specific configurations providing a strong 
field enhancement such as, for example, resonant micro- 
cavities [6-8] or photonic crystals [9-12]. Analogously, 
the field enhancement occurring on a rough metal sur- 
face due to the excitation of surface plasmon polari- 
tons is responsible for a substantial enhancement of the 
surface SHG [13]. Conceiving subwavelength coherent 
light sources is a fundamental issue of modern nanopho- 
tonics [13, [lH so that achieving SHG from nanostruc- 
tures is an important target. At the nanometer length 
scale, the small interaction distances generally entail an 
highly inefficient SHG so that peculiar mechanisms for 
locally enhancing the electromagnetic field of the fun- 
damental wave have to be harnessed. As an example, 
the plasmonic field enhancement [16] occurring in vicin- 
ity of the subwavelength apertures in a metallic film (ac- 
companying the extraordinary linear optical transmission 
[l7]) is responsible for an efficient SHG (l8l42lj when the 
holes are filled by a quadratic nonlinear material. An- 
other strong plasmonic field enhancement occurs within 
spherical nanocavities with dielectric core and plasmonic 
nanoshell [22] so that, if the core is filled with a noncen- 
trosymmetric nonlinear medium, a lager enhancement of 
the SHG is observed 

In this Letter we theoretically show that a subwave- 
length thick slab consisting of a quadratic nonlinear 
medium whose linear permittivities are very small and 
have different signs is able to provide highly efficient 
SHG. Properties and applications of epsilon-near-zero 
materials have recently been proposed both in the lin- 
ear [23-^28'] and in the nonlinear regime [29'-^33] and in- 
definite materials [34-37] have attracted a good deal of 
attention as well. Here we point out that, the small 
value of the dielectric permittivities produces a strong 




FIG. 1: (Color online) Slab and fields geometry for the SHG 
process. 



non-plasmonic enhancement of the longitudinal field per- 
pendicular to the slab surface. The field enhancement is 
simply a consequence of the fact that, in conditions of 
linear complete slab transparency, the longitudinal fun- 
damental field within the slab across the vacuum-slab 
interface coincides with the longitudinal field of the inci- 
dent fundamental wave divided by the very small dielec- 
tric permittivity. We predict that the enhanced longitu- 
dinal field is responsible for a SHG process with efficiency 
up to 40%, for impinging optical intensities of a few of 
MW/cm^. 

The geometry of the considered SHG setup is sketched 
in Fig.l and, labelling hereafter the fundamental and 
second-harmonic quantities with the superscripts (1) and 
(2), we choose A^^) = 827 nm and A^^) = A(^V2 = 
413.5nm for the two wavelengths. The slab of thickness L 
consists of quadratic nonlinear anisotropic medium whose 
relative dielectric tensor is e = diag (e^^, e^, e^) (i.e. its 
principal axes coincide with cartesian axes) and whose 
second-order nonlinear optical tensor has the only non- 
vanishing components di4 = ^25 = <^36 = x/2. Here we 



choose ei^^ 



-0.01, ei^^ = 0.001 and e\f^ = 2.86 for the 



only permittivities playing a role in our discussion and 
X = 1.3pm/V for the nonlinear susceptibility (see below 
for a possible composite medium exhibiting these prop- 
erties). The slab is illuminated by a plane fundamental 
wave (FW) with incidence angle and with its electric 
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FIG. 2: (Color online) Logarithmic plot of the linear slab 
transmissivity T (panel a) and of the linear longitudinal field 
enhancement factor F (panel (b)) for various incidence angles 
and slab thicknesses L. The solid lines along the bottom 
planes and marked with n — 1,2,3 are the curves of linear 
complete slab transparency (i.e. K'^^^ L — rnr). 



field polarized in the plane of incidence (p-polarization) 
so that the reflected (r) and transmitted (t) plane FWs 
are p-polarized as well. The slab also scatters reflected 
(r) and a transmitted (t) plane second harmonic waves 
(SHW) whose electric fields are polarized perpendicular 
to the plane of incidence (s-polarization), as a conse- 
quence of the chosen anisotropic quadratic nonlinearity 
(see Fig.l for the definitions of the field amplitudes and 
wave vectors). 

As a prelude to the analysis of the SHG process, it 
is worth discussing the linear slab behavior observable 
when the optical intensity of the incident FW is very 
small. The linear slab transmissivity is 
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where K^^\0) = k^^^ e^^\l - sin^ is the 

wave vector of the radiation excited within the slab, 
A:(i) = 27r/A(i) and F = k^^h^^\os / (2K^^^) + 

K^'^y (2k^^h^x^ cos l9). In Fig.2(a) we plot the logarithm 
of T as a function of both the incidence angle and the 
slab thickness L from which we note that the slab is prac- 
tically always opaque to the considered radiation except 
for very specific incidence angles (dependent on the slab 
thickness) for which it is completely transparent. From 
Eq.([T]) it is evident that total transmission T = 1 occurs 
if K^^^L = nir and the ensuing curves of complete trans- 
parency are reported with solid lines along the bottom 
plane of Fig.2(a). Note that if K^^^L is slightly differ- 
ent from nTT, the denominator in Eq.([T]) is dominated by 



the term containing F (since \F\ 



.(1)1 



> 1) so that 



the overall transmissivity turns out to be proportional 
to lei^^l thus explaining the sharpness of the ridges of T 



of Fig. 2 (a) corresponding to complete transparency. In 
Fig. 2(b) we plot the logarithm of the longitudinal field 

enhancement factor F = le^ • E(^)(0+)|/|e;2 • E-"^^| defined 
as the ratio between the modulus of the field perpen- 
dicular to the slab within the medium just across the 
interface z = (see Fig.l) and the modulus of longitu- 
dinal field of the sole incident FW. Note that, in corre- 
spondence of the curves of complete slab transparency, 
the field enhancement factor is F = 1000 = l/ei^\ 
This is easily understood from one of the field match- 
ing condition at the interface 2: = (continuity of nor- 
mal component of the electric displacement field, i.e. 



(1) 



E^^^) = e^'^[e, • E(i)(0+)]) since, when T = 1 



.(1) 



no reflected wave is generated (i.e. Er"*^^ = 0), so that 

•E(i)(0+) = [e^ -Ef V^i^^. We conclude that the con- 
sidered slab provides a powerful mechanism for attaining 
a strong enhancement (amounting to 1/ei^^) of the lon- 
gitudinal field. Note that, for the field enhancement to 
occur it is essential that i^(^) is real and this motivates 
our choice of the permittivities' signs. 

In order to discuss the SHG process we note that the 
incident p-polarized FW (see Fig.l) produces, within 
the slab, the transverse magnetic (TM) fundamental 
field E(i) = e^^^'^^^i"^[^i'\z)e^ + ^^'^(^)e,], H(i) = 
e'^^^^^^^^^^[Hy^\z)ey] which, due to the chosen quadratic 
nonlinearity, couples to the transverse electric (TE) 
second-harmonic field E(2) = e^2^^'^^^^"^[£;f ^(z)e^]. 



H(2) = ei2fe^ ^ ^^^^^[Hr {z)ea: ^ Hr {z)e^] (here we have 
assumed e"*'^^ and e"*^'^^ as time factors for the FW 
and the SHW fields, respectively). The second-harmonic 
TE field within the slab matches with the s-polarized 
reflected and transmitted plane SHWs scattered by the 
slab [38]. The dielectric slab response is described, in the 
considered situation, by the electric displacement fields 

for the FW and 
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for the SHW. After setting = 1/2 Ve^T^olEf ^ p = 
7.9MW/cm^ for the optical intensity of the incident FW, 
we have numerically solved, for various incidence angles 
and slab thicknesses L, the coupled sets of Maxwell 
equations (within the slab < z < L, where Eqs.(j2j) and 
© hold) for the FW and the SHW, along with match- 
ing conditions (at z = and z = L) describing the slab 
illumination by the sole incident FW. Consequently, we 
have evaluated the slab efficiencies r]r = |Er^^ p/lE^"*^^ p 



and r]t 



ie: 



(2),2/|e(1)|2 



of converting the incident FW 



into reflected (for 2: < 0) and transmitted (for z > L) 
SHWs and their logarithm plots are reported in Fig. 3 (a) 
and 3(b), respectively. From Fig. 3, it is particularly ev- 
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FIG. 3: (Color online). Logarithmic plots of the slab efficien- 
cies ?7r (panel (a)) and rit (panel (b)) of converting the inci- 
dent FW into reflected and transmitted SHWs, respectively. 
The SHG efficiencies are evaluated for the optical intensity 
/^^) = 7.9MW/cm2 of the incident FW. The solid lines along 
the bottom planes and marked with n — 1, 2, 3 are the curves 
of linear slab complete transparency (i.e. K'^^^ L = nvr). 




FIG. 4: (Color online) Values of the SHG efficiency rft (ex- 
tracted from the numerical simulation of Fig. 3) for incident 
angles and slab thicknesses corresponding to the linear slab 
complete transparency for the FW. 



ident that, for the chosen optical intensity of the inci- 
dent FW, SHG uniquely occurs for specific incident an- 
gles and slab thicknesses practically coinciding with those 
providing linear total transmission of the FW (see the 
sharp ridges of rft and rfr approximately located upon 
the curves K^'^^L = nir). The physical origin of such 
phenomenology is easily grasped by noting that SHG is 
here due to the source term proportional to the com- 
ponents Ex^"^ and E^^"^ of Eq.(|3]), so that, due to the 
small value of the nonlinear susceptibility the SHW 
field scattered by the slab is generally very weak for the 
chosen impinging optical intensity I^^^ In these condi- 
tions, the nonlinear terms in Eqs.(|2]) can be generally 
neglected so that the FW does not experience the ef- 
fect of the quadratic nonlinearity and it consequently ex- 
hibits the linear behavior of Fig. 2. Such an unattrac- 
tive scenario is fundamentally altered when the condi- 



tions of linear complete slab transparency for the FW 
are met, since the above described enhancement mecha- 
nism of the longitudinal field E^^"^ entails an enhancement 
of the source term x^i^"^ E^^"^ so marked to yield efficient 
SHG (see the sharp ridges of Figs. 3 (a) and 3(b)). In 
Fig. 4, we plot the values of the above evaluated SHG 
efficiency r]t exclusively along the curves of linear total 
transmission of the FW. From Fig. 4 we conclude that 
the proposed mechanism is very efficient, even at the 
chosen optical intensity /-^^^ = 7.9 MW/cm^ since it is 
characterized by values of the efficiency r]t up to 40% 
and it takes place in slabs with thickness spanning the 
subwavelength range < L < 500 nm. For the sake 
of completeness, note that the proposed SHG scheme 
does not rely on the phase-matching condition since 

the slab SHW wavevector is K^'^^ = 2k^'^^ ^ ef^ - sin^ l9 
(in the undepleted pump regime) so that the condition 
2i^(^) = i^(^) yields the phase matching angle 6pm = 

asin^ei^^(e^^^ - ef^)/{e;t^ - ei^^) 30.35°, an angle not 
playing any special role in the above described process 
(see Figs. 3 and 4). 

The above analysis has been carried out for a specific 
value of In order to discuss the role played by the 
FW intensity on the considered SHG process, we have 
performed further numerical simulations in a L = 400 nm 
thick slab to evaluate the SHG efficiency as a function of 
the incidence angle and the optical intensity The 
logarithm of the obtained rjt is plotted in Fig. 5 (a) (where 

the axis is logarithmic) and we have considered a 
restricted angular range (37° < < 45°) surrounding 
the angle = 40.9° which is one the angles (the one be- 
longing to the n = 2 branch) at which the 400 nm thick 
slab is linearly transparent to the FW (see Fig. 2). Note 
that, for angles not close to = 40.9°, the efficiency rjt 
is a growing function of whereas, at 6> = 40.9°, the 
efficiency grows for intensities up to 50 MW/cm^, it un- 
dergoes a saturation and eventually decreases at higher 
intensities. For clarity purposes, in Fig. 5(b) we plot the 

efficiency rft as a function of ^ for = 40.9°. Such a be- 
havior of r]t is evidently a consequence of the fact that, at 
low optical intensity, the undepleted pump regime rules 
the phenomenology and provides the growing efficiency 
rjt whereas, at higher intensities, the full nonlinear regime 
comes into play where the FW and the SHW experience 
their mutual nonlinear interaction. It is however remark- 
able that, in the situation considered in Fig. 5, the max- 
imum efficiency rjt = 0.316 is achieved at the intensity 
/^(i) ^ 55 MW/cm^ (see Fig.5(b)) which is considered to 
be a small value in a nonlinear optical setup. 

In order to prove the feasibility of the above described 
SHG process, we discuss here the numerical design of 
a nano-structured composite characterized by the above 
unconventional features. The slab (see. Fig.l) is man- 
ufactured by alternating, along the axis, two kind of 
layers (labelled with a and b) of thicknesses da = 11.84nm 
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FIG. 5: (Color online) Efficiency rft of a 400 nm thick slab of 
converting the incident FW into the transmitted SHW. (a) 
Logarithmic plot of rft as a function of the incidence angle 
and the optical intensity of the incident FW. (b) Linear 
plot of r]t dit = 40.9° (detailed from the shadowed section of 
panel (a)). 



and db = 88.16 nm, respectively. We choose the lay- 
ers of kind a to be filled with alpha- iodic acid (a-HIOs) 
crystals with principal axis along the cartesian axis so 
that the permittivities relevant for our purposes are 



^ax — 3.7588, CaJ 



3.8517, e^aj = 3.5871 and the non- 



vanishing components of the second-order nonlinear op- 
tical tensor are di4 — ^25 = <^36 = Xa/2 = 7 pm/V 
[39]. For the layers of kind h we choose a polymer (e.g. 
poly (methyl methacrylate), PMMA) with permittivities 

ty^^ = 2.2022 and ^ = 2.2575 hosting silver-coated 
spherical PbSe quantum dots (QDs) with inner radius 
Ri = 5 nm and outer radius R2 = 6.38 nm. Exploiting 
the electrostatic approach of Ref. [40] , the structured par- 
ticle can be described by the equivalent permittivity es = 
eAg[eQD{l + 2p) + 2eAg{l - p)]/[eQD{l - p) + 2eAg{2 + p)] 
where p = caq = e 00— ujp / {uj'^ -{-iTou) is the Drude- 

type silver permittivity (eoo = 4.56, ujp = 1.38 • 10^^ s~^ 



and r = 0.1 • 10^^ s-i [41]) and cqd = e& + ^^o/(^ - 
uJq -\- 12^00) is the dielectric permittivity of the considered 
semiconductor QD (65 = 12.8, ujq = 2.2793 • 10^^ s-\ 
7 = 1.2788 • 10^2 s-i [42] and A = 9.6- 10"^). Assuming 
Fqd =0.1 for the volume filling fraction of the QSs em- 
bedded in the polymer and using the standard Maxwell 
Garnett mixing rule for obtaining the effective dielectric 
permittivity, we obtain, for the layers of kind 6, e^^^ = 

-0.5160 + 0.0002Z and ef^ = 2.7819 + 0.0127z. Since the 
period da-\-db = 100 nm is considerably smaller than the 
two wavelengths, FW and SHW propagating through the 
slab will experience the effect of a homogeneous medium 

characterized by the permittivities e^^^ = fa^ax-\-fb^\^^ = 
(1) _ . .(1) ^ ^^,(1) 



-0.0099 + O.OOOlz, e 



O.OOOH, e^y^ 



.(2) 



= fa 

h/e 



0.0011 



(2)i 
h 



= 2.8579 + 0.01184Z 

(where fa = da /{da + 4) and fb = 4/(<^a + 4) are the 
layers volume filling fractions), whereas the effective sus- 

ceptibilityjs x = fa{ef^ [eay)Xa = 1.3206+0.0054zpm/V 
(see Ref.[43| for the derivation of these relations). Even 
though theoretical (note that we have neglected the un- 
certainties of all the used quantities), the proposed ap- 
proximate design reveals that, by suitably acting on free 
parameters (kind of media, QDs volume filling fractions, 
inversion population factor A, etc.), a medium for observ- 
ing the proposed efficient SHG can actually be conceived. 

In conclusion, we have identified a novel mechanism 
able to provide a strong field enhancement in anisotropic 
optical metamaterials having very small cartesian per- 
mittivities of different signs. We have showed that, if 
the slab has a nonlinear quadratic response, the field 
enhancement is responsible for an efficient SHG process 
even if the slab thickness is smaller than the wavelength 
of the SHW. It is worth stressing that the proposed field 
enhancement mechanism can easily be observed also in 
setups different from the considered slab. Our findings 
offers an hitherto unexplored way for achieving efficient 
SHG from nanostructures and therefore we believe they 
could pave the way for conceiving a novel generations of 
nanometric-sized coherent light sources. 
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